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Abstract 

We consider a power-controlled wireless network with an established network topology in which 
the communication links (transmitter-receiver pairs) are corrupted by the co-channel interference and 
background noise. We have fairly general power constraints since the vector of transmit powers is 
confined to belong to an arbitrary convex polytope. The interference is completely determined by a so- 
called gain matrix. Assuming irreducibility of this gain matrix, we provide an elegant characterization of 
the max-min SIR-balanced power allocation under such general power constraints. This characterization 
gives rise to two types of algorithms for computing the max-min SIR-balanced power allocation. One 
of the algorithms is a utility-based power control algorithm to maximize a weighted sum of the utilities 
of the link SIRs. Our results show how to choose the weight vector and utility function so that the 
utility-based solution is equal to the solution of the max-min SIR-balancing problem. The algorithm is 
not amenable to distributed implementation as the weights are global variables. In order to mitigate the 
problem of computing the weight vector in distributed wireless networks, we point out a saddle point 
characterization of the Perron root of some extended gain matrices and discuss how this characterization 
can be used in the design of algorithms in which each link iteratively updates its weight vector in parallel 
to the power control recursion. Finally, the paper provides a basis for the development of distributed 
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power control and beamforming algorithms to find a global solution of the max-min SIR-balancing 
problem. 

Index Terms 

Max-min SIR-balancing, Max-min fairness, Power control, Wireless networks, Utility maximization, 
Interference management, Distributed algorithms 

I. Introduction 

Wireless channel is error-prone and highly unreliable being subject to several impairment 
factors that are of transient nature, such as those caused by co-channel interference or multipaths. 
Excessive interference can significantly deteriorate the network performance and waste scarce 
communication resources. For this reason, strategies for resource allocation and interference 
management are usually necessary in wireless networks to provide acceptable QoS levels to the 
users. 

There are different mechanisms for resource allocation and interference management. Power 
control may play a central role in distributed wireless mesh networks, where, due to the lack 
of a central network controller, link scheduling strategies are notoriously difficult to implement. 
Thus, a reasonable approach is to avoid only strong interference from neighboring links, and 
then use an appropriate power control policy to manage the remaining interference in a network. 
In this paper, we focus on the power control problem, which addresses the issue of coordinating 
transmit powers of links such that the worst signal-to-interference ratio (SIR) balanced against 
some SIR targets attains its maximum. This so-called max-min SIR-balancing problem is a 
widely studied resource allocation problem for wireless networks (see, for instance, fl], ED, H, 
El, 0, & 0, [H, @, OH, OH, EH as well as EES Sections 3.1 and 5.6], El, 031 Section 
5.6]). A key feature of this strategy is that any given SIR (signal-to-interference ratio) targets 
are feasible if and only if they are satisfied under a max-min SIR-balanced power allocation. 
Moreover, the notion of max-min SIR-balancing is closely related to max-min fairness, the 
most common notion of fairness. Note that since we focus on static wireless networks, transmit 
powers are to be periodically adjusted to changing channel and network conditions (dynamic 
power control). This in turn presumes a relatively low up to moderate network dynamics. In 
contrast, in highly dynamic wireless networks, one should consider resource allocation schemes 
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for stochastic wireless networks IfTBH, |fT7l, IfTSl, lfT9l. GUI. IZD. E21. 11251. 

In the noiseless case, which is widely considered in the literature [QQ|, ED, J51, 0, 01, ll24l . 
Il25l . [|26l . (an overview can be found in IfTTTl . Il27l ) and where power constraints play no role 
in the analysis, it is widely known IfTTTl . |[T3l Sections 3.1] that any positive eigenvector of the 
(irreducible) gain matrix scaled by a diagonal matrix of given SIR targets is a solution to the 
max-min SIR-balancing problem. In [12J, the problem was solved for a "noisy" downlink channel 
constrained on total power. The sum constraint on transmit powers was captured by an additional 
equation so that the optimal solution is characterized in terms of some unique eigenvector of a 
certain irreducible gain matrix of higher dimension (see also |[T3l pp. 111-113]). 

Assuming an irreducible gain matrix, Sections ITITI - ITVl extend these results to any convex 
poly tope as the constraint set to model constraints on transmit powers of the links. In addition 
to the analysis in a higher dimension (as in 021), we also obtain an elegant characterization 
of the max-min SIR-balanced power vector from an eigenvalue problem of the same dimension 
as the original problem. These results were inspired by ll28l . where the authors used different 
tools to characterize rate region in interference channel with constrained power (see also the 
acknowledgments after the conclusions at the end of this paper). 

In Section |Vj we use the results of Sections ITITMlVl to establish a connection between the 
max-min SIR-balancing power control problem and the utility-based power control problem. 
Such a connection is known in the noiseless case lfT5l Section 5.9] and constitutes the starting 
point for the analysis in lfl4l . More precisely, we show how to choose the weight vector and 
utility function so that solving the problem of maximizing a weighted sum of the utilities of 
the SIRs leads to the max-min SIR-balancing solution. This result was used in ll29l to solve the 
max-min SIR-balancing problem over the joint space of admissible power vectors and receive 
beamformers. Thus, the results of this paper provide key tools to generalize the results of lfl4l 
to noisy channels under general power constraints and to a larger class of utility functions. 

An advantage of the utility-based approach is that there exist distributed power control schemes 
to compute the max-min SIR power vector, provided that each link knows how to select its weight 
ll30l . OTI . The problem is however that a desired weight vector is determined by positive left 
and right eigenvectors of some nonnegative matrix, so that the links cannot choose their weights 
independently. Thus, as neither the eigenvectors nor the corresponding matrices are a priori 
known at any node, the presented approach for computing the max-min SIR power allocation is 
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still not amenable to implementation in decentralized wireless networks. 

A basic idea to overcome or at least to alleviate this problem is to let each link iteratively 
update its weight vector in parallel to the power control recursion. In Section IV-B[ we point 
out that a saddle point characterization of the Perron roots of some nonnegative matrices. This 
characterization provides a basis for efficient saddle point algorithms converging to a max- 
min SIR-balanced power vector. Basically, the idea is redolent of primal-dual algorithms that 
employ some optimization methods to minimize the Lagrangian over primal variables and to 
simultaneously maximize it over dual variables. 

Before starting with the analysis, we introduce definitions, notation and state the max-min 
SIR-balancing problem. 

II. Definitions and problem statement 

We consider a wireless network with an established network topology, in which all links share 
a common wireless spectrum. Let K > 2 users (logical links) compete for access to the wireless 
links and let X = {1, ... , K} denote the set of all users. The transmit powers p k , k £ X, of the 
users are collected in the vector p = (pi, . . . ,Pk) > 0, which is referred to as power vector or 
power allocation. The transmit powers are subject to power constraints so that p G P wherd_ 

P = { p G : Cp < p, C G {0, 1} N * K } C R K (1) 

for some given p = (Pi, . . . , P N ) > and C with at least one 1 in each column so that P is a 
compact set. Throughout the paper, we use N = {1, . . . , N} where N is the number of power 
constraints. The main figure of merit is the SIR at the output of each receiver given by 
(A.l) SIRfc(p) = pk/h(p), k E X, where the interference function I k is h(p) = (Vp + z) k = 

sr^K 

1^1=1 V k,lPl + Zk- 

V := (fjtj) € R+ xK is the gain matrix, v k ,i = V^i/V^k if Z ^ A; and if I — k where Vf-j > 
with Vk t k > is the attenuation of the power from transmitter / to receiver k. The A;th entry of 
z := {z\, . . . , z K ) is z k = Ofc/Vfc^, where o\ > is the noise variance at the receiver output. 

Let 7i, . . . , 7# > be the SIR targets and let T = diag(7i, . . . , 7^). We say that the SIR 
targets are feasible if there exists a power vector p G P (called a valid power vector) such that 
SIR fe (p) > 7fc > 0. 



'R+,R++ are nonnegative and positive reals, respectively. 



5 



Definition 1: Given any T, p is said to be a max-min SIR-balanced power vector if 

p := argmaxmin (SIR fe (p)/7 fc ) . (2) 

It is important to notice that the maximum in © exists as minfc e :x:(SiRfc(p) /jk) is continuous 
on the compact set P. Thus, 7^ are not necessarily met under optimal power control ©. For 
this reason, 7^ can also be interpreted as a desired SIR value of link k. A trivial but important 



observation is that p > 0, allowing us to focus on P + = P n 



III. Some Preliminary Observations 

In general, p of Definition Q] is not unique. For general power constraints, the uniqueness is 
ensured if V is irreducible 11321 . Il33l since then the links are mutually dependent through the 
interference. In order to see this, notice that the problem © is equivalent to finding the largest 
positive threshold t such that t < SIR/^p)/^ for all k G % and p G P. The constraints can 
be equivalently written in a matrix form as Tz < (|l — TV)p and p G P. So, as p must be a 
positive vector, lfl"5l Theorem A.51] implies that the threshold t must satisfy 

p(rv)<i/t. O) 

Now, one particular solution to the max-min SIR-balancing problem © is p' given by 

p' = (i/t'i-rv) _1 rz,p' g p (4) 

where 

t' = arg max t > 1 s.t. (l/tl — TV)p = Tz, p G P . (5) 

Note that p' is a max-min SIR-balanced power vector such that SlRk{P')/lk = SIRj(p')/7z for 
each k, I G %. This immediately follows from ([4]) when it is written as a system of K SIR 
equations. By ©, ©, © and lfl"5l Theorem A. 5 11, p' > exists and is a unique power vector 
corresponding to a point in the feasible SIR regioro that is farthest from the origin in a direction 
of the unit vector t/ ||t II 1- 

The above considerations are illustrated in Fig. \T\ The plots depict two examples of feasible 
SIR regions in a system with two users. In both cases, the point (SIR(p')i, SIR(p') 2 ), with p' 

2 The feasible SIR region is the subset of R^: of all SIR levels that can be achieved by means of power control. F defined by 
10 becomes the feasible SIR region if <j>(x) = x and jk = 1, k G DC. 
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Fig. 1. The feasible SIR region under individual power constraints and two different gain matrices V > 0. The following 
notation is used: 7^ = SIRfc(p) and 7^ = SIRfc(p') where p and p' are defined by l|2} and {4j, respectively. Left: V is chosen 
so that SIR2(p) = P2/Z1, in which case p is not unique. Right: V is irreducible, in which case p is unique and equal to {4j. 
The point (7", 72) corresponds to the max-min fair power allocation. 



given by © is the point at the intersection of the boundary of the feasible SIR region with 
the line defined by the 7 vector and is max-min SIR-balanced. It is however not the unique 
solution if SIRs of some users can be increased without affecting the minimum. The examples 
of configurations in which the power allocation p' given by © is not and is a unique solution 
of the max-min SIR-balancing problem are presented in the left and in the right subplot of Fig. 
[TJ respectively. 

Now let us assume throughout the paper that 

(A.2) : M ++ — > Q C R is continuously differentiable and strictly increasing function. 

By strict monotonicity, one has 0(min fcg3 <;SIR fc (p)/7 fc ) = min fce3C 0(SIR fc (p)/7 fc ) for every 

p > 0. Thus, as p G P+, 

p = argmaxmin0(SIR fc (p)/7 fc ) (6) 
peP+ kex 

where p is a max-min SIR-balanced power vector defined by ©. With © in hand, we can 
prove a sufficient condition for the uniqueness of p under general power constraints. To this 
end, given 0, we define the set F C as 

F = {q G Q K : q k = 0(SIR fc (p)/ 7fc ) , k G X, p G P+} (7) 

Note that F can be interpreted as a feasible QoS region where the QoS value for link k is 
0(SIR fc (p)/7 fc ). By |(A.2)| and lfT51 Sect. 5.3], the following can be said about F. 

Observation 1: There is a bijective continuous ma^\ from F onto P + (see also (1231)). If 



3 The reader should bear in mind that the existence of such a bijective map allows us to prove some results in F. 
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(A.3) the inverse function g(x) := 1 (x) is log-convex |fT5l , 

then F is downward comprehensive, connected and convex. Moreover, SIRfc(p') with p' defined 
by © corresponds to a point on the boundary of F where q x = . . . — qn with q k defined by ©. 

Note that the boundary of F (denoted by <9F) is the set of all points of F such that, if 
p is the corresponding power vector in ©, then Cp < p holds with at least one equality. 



Widely known examples of functions satisfying |(A.2)| and |(A.3)| are x t— > \og(x),x > 0, and 



x \— > — l/x n ,n > l,x > 0. Now let us state the following auxiliary result. 



Lemma 1: Suppose that |(A.2)| and |(A.3)| hold. Let V be irreducible. Then, q G <9F if and 
only if there exists w > such that q maximizes x i— > w T x over F. 

Proof: By convexity and downward comprehensivity of F, every boundary point of this set 
is a maximal point and it maximizes x i— > w T x over F for some w > [|34l pp. 54-58]. Since 
V is irreducible, it follows from ll35l Theorem 4.3] (see also ll35l Corollary 4.3]) that w > 0. ■ 

Now we can easily observe the following. 

Observation 2: If V > is irreducible, then p is unique and equal to p' defined by ([4]). 
Proof: The proof is deferred to Appendix IVIII-AI ■ 

We complete this section by pointing out a connection between the max-min SIR-balancing 
problem considered in this paper and the notion of max-min fairness. A vector of balanced 
SIRs with entries SIR fc (p")/7 fc , k G X is max-min fair if any SIR fc (p")/7 fc cannot be increased 
without decreasing some $IRi{p")/yi, I ^ k which is smaller than or equal to SIRk(p") / Ik', the 
vector p" is then a max-min fair power allocation ll36l . |[37ll . A max-min fair power allocation 
is therefore also max-min SIR-balanced (provided that the feasible SIR region is downward 
comprehensive); the converse is in general not true. However, if the max-min SIR-balancing 
problem has a unique solution given by ©, this solution is also max-min fair and there are 
no other max-min fair power allocations. These relations can be observed in Fig. [T] where both 
max-min SIR-balanced and max-min fair points are indicated. 

IV. Characterization of Max-Min SIR-balancing 

In this section, we characterize p G P + defined by © under the assumption that V is 
irreducible. We point out possible extensions to reducible matrices at the end of this section. 
We assume that P C IR+ is a convex polytope given by Q). So, throughout this section, 
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max nGJ4 g n (p) < 1 where 

g n (p) := 1/P n c£p, tiGN (8) 

and c n G {0, 1} A is a (column) vector equal to the nth row of the matrix C. Now, using ([8]), 
the max-min SIR power vector p defined by © can be written as 

p = argmaxmin(SIR A; (p)/7 fc ) s.t. max# n (p)<l. (9) 

where 7^ > 0, k G X, is arbitrary but fixed. 

Lemma 2: Let p be any power vector that solves ©. Then, the following holds 

(i) max n6 ^5' n (p) = 1. 

(ii) If V > is irreducible, then p is unique and 

V fceaC 7fc/SIR fe (p) = /5, /3>0. (10) 

Proof: The proof can be found in Appendix IVIII-B1 ■ 
Because p maximizes min fcg3 <;(SIRfc(p)/7fc) over P, it follows from (flOl) that 1/(3 > is the 
corresponding maximum. It must be emphasized that (flOl) is not true for general nonnegative 
matrices V > 0. In the lemma, we require that the gain matrix be irreducible, which is sufficient 
for (flOl) to hold but not necessary. The irreducibility property ensures that, regardless of the 
choice of P, there is no subnetwork being completely decoupled from the rest of the network. 
To be more precise, if V is irreducible, then the network is entirely coupled by interference 
so that the type of power constraints is irrelevant for this issue (see also the remark at the end 
of this section). Unless otherwise stated, it is assumed in the remainder of this section that 
V > is an arbitrary irreducible matrix. Due to (ii) of Lemma [2l this implies that the max-min 
SIR-balanced power vector is unique. 
Let us define 

N (p) := {^0 e N : n = argmax# n (p) = l} (11) 

which includes the indices of those nodes for which the power constraints are active under the 
power vector p. By (i) of Lemma [2l the cardinality of N (p) must be larger than or equal to 1. 
In what follows, let (3 > be the constant in part (ii) of the lemma. This together with part (i) 
implies that 

/?p = rvp + r z <?n(p) = i- (12) 
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Putting the first equation into the second one yields a set of K + 1 equations that, if written in a 
matrix form, show that if p solves the max-min SIR-balancing problem, then there is a constant 
P > such that 

pp = A^p, p > 0, p G (13) 

for each n G N (p) where p = (p, 1) is the extended power vector and the nonnegative matrix 
A W e Rf +1 )x^ +1 ) i s defined to be 

, , / rv r z \ 

A (n) = , n£N. (14) 

U c " rv ^ Tz ) 

Alternatively, we can write (fT2)) as (3p = TVp + Tz ■ g n (p), from which we obtain, for each 

n G 3ST (p), 

/?p = B^p, /3>0,pGM^+ 1 (15) 
where B (n) eRf K is defined to be (for each n G X} 

b (b) := rv + -Lrzc^ = r(v + i-zcj) = rv (n) ( i 6) 

± n i n 

and V := V + -p-zc^,n G N. So, given p defined by ©, the equations <fT3l and (IT5T ) 
hold for each n G 3sT (p). In other words, the solution of © in a network entirely coupled by 
interference must satisfy (fT3l) and (fT5l) for each node n G N whose power constraints are active 
at the maximum. This is summarized in the following lemma. 

Lemma 3: If V > is irreducible and p solves the max-min SIR-balancing problem ©, then 
p satisfies both (fT3l) and (fT5l) for some (3 > and each n G N (p). 

Note that the lemma is an immediate consequence of parts (i) and (ii) of Lemma [2l from 
which (fT3l) and (fl~5T) follow for an arbitrary n G N (p). Now we are in a position to prove the 
following result. 

Lemma 4: Suppose that V > is irreducible. Then, for any constants c\ > and C2 > 0, the 
following holds. 

(i) For each n G N, there is exactly one positive vector p = p( n ) G with g n (p) = ci 
satisfying /?( n ) p = B^p for some (3^ > 0. Moreover, (3^ is a simple eigenvalue of 
B (n) and = p(B (n) ). 
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(ii) For each n £ INT, there is exactly one positive vector p = p^ n > £ with p^+i = c 2 

satisfying /3( n ) p = A^p for some > 0. Moreover, (3^ is a simple eigenvalue of 

A (n) and pW = p(A (n) ). 

Proof: The reader can find the proof in Appendix IVM-CI ■ 
The lemma says that, for each n £ N, the matrix equation /?( n ) p = B^p with f3^ > 
and p £ IR+ is satisfied if and only if p is a positive right eigenvector of B^ associated with 
p( n ) = p(B^). Furthermore, if g n (p) = L then p is unique. Similarly, (5^ p = A^p with 
pi n ) > o and p £ IR+ +1 is satisfied if and only if p is a positive right eigenvector of A^ 
associated with (3^ = p(A ( - n ' ) ) and there is exactly one such an eigenvector whose last entry 
is equal to one. Furthermore, for each n £ X, (3 {n) = p(A (n) ) = p(B (n) ). This is because 
if p(B (ri) )p = B (n) p holds for some n £ N, then we must have p(B (n) )p = A (n) p where 
p = (p, 1) £ Thus, by Lemma HI we must have 

p(A (n) ) =p(B (n) ), neX. (17) 

Note that a solution to the max-min SIR-balancing problem is not necessarily obtained for each 
n £ N since, in the optimum, some power constraints may be inactive. Indeed, in general, the 
set Nq(p) = X \ No(p) is not an empty set where p defined by © is unique due to (ii) of 
Lemma [2] 

Now we combine Lemmas [3] and @] to obtain the following. 

Theorem 1: Let [3 be given by (flOl) . If V > is irreducible, then the following statements 
are equivalent. 

(i) p £ P + solves the max-min SIR-balancing problem ©. 

(ii) For each n £ N (p), p is a unique positive right eigenvector of B^ associated with 
P = p(B {n) ) > such that g n (p) = 1. 

(iii) For each n £ N (p), p is a unique positive right eigenvector of A^ associated with 
/3 = p(A (n) ) > such that p K+1 = 1. 

Proof: The proof can be found in Appendix IVIII-Dl ■ 
Theorem Q] implies that if V is irreducible, then p > is the (positive) right eigenvector of 
B (n) associated with p(B (n) ) £ <r(B (n) ) for each n £ N (p). Alternatively, p can be obtained 
from p = (p, 1), which is the positive right eigenvector of A^ associated with p(A^) for each 
n £ 3Sfo(p). The problem is, however, that N (p) is not known as this set is determined by the 
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solution to the max-min SIR-balancing problem, and hence its determination is itself a part of 
the problem. As the SIR targets are feasible if and only if they are met under p, the following 
characterization of the set N (p) immediately follows from [|28l and (TT71) . 
Theorem 2 ([28]): Suppose that V > is irreducible. Then, 

N (p) = {n G N : n = argmax neN p(BW)} (18) 
Moreover, the feasible SIR region F 7 is characterized as 

F 7 ={ 7 6 : maxp(B (n) ) = maxp(A (n) ) < 1, 

r = diag( 7 i,...,7A')}. (19) 

The characterization of the feasible SIR region in (TT9l) can be deduced directly from ll28l as 
the authors show that min(SIRfc(p)/7fc) < min(l/p(B n )), without characterizing, however, the 
corresponding power allocation vector p. 

Theorems Q] and [2] directly lead to the following procedure for computing the max-min SIR 
power vector p given by ©. 

Input: T = diag(7i, . . .,7^). 
Output: peP + cP 

l: Find an arbitrary index n G N such that n = argmax ne:N - p(B^) where B 1 -™^ is given by 
©. 

2: Let p be given by p(B^ n °^)p = B^p and normalized such that c^ p = P no . 



Remark 1 (Remark on reducible matrices): We point out that all the statements in this paper 
hold if B < - n ' ) is irreducible for each n G N, which may be satisfied even if V is reducible. This is 
for instance true in the case of a sum power constraint (C = (1, . . . , 1)) where B = B^, 74 = 
{1}, is a positive matrix. Moreover, the statement of Theorem Q] can be shown to be true even 
if B 1 -™) is irreducible only for n G N (p). 

V. Applications 

In this section, we discuss two other applications of the results. In doing so, V is assumed 
to be irreducible. Under this assumption, the feasible QoS region given by CD) can be shown to 
be strictly convex ll35l Corrolary 4.3]. 
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A. Computation via utility-based power control 

In this section, we show that p can be obtained by maximizing the following aggregate utility 



function 



F(p,w) := V W>(SIR fc (p)/ 7fe ) 
provided that the weight vector w = (w±, . . . , wk) > is chosen appropriately and 



(20) 

-> Q 



satisfies |(A.2)| and |(A.3)| To this end, define 

p*(w) := argmax peP+ F(p, w) . (21) 
Although P + is not compact, it can be shown |fT5l that the maximum exists if|(A.2)|and|(A.3) 



are 



fulfilled. Furthermore, it is obvious that in the maximum, at least one power constraint is active, 
that is, Cp* < p holds at least with one equality. Thus, we have q*(w) = (<Zi(w), . . . , q* K {vr)) G 
dF for 

q* k (w) = 0(SIR fc (p*(w))/ Tfc ) • (22) 

In words, p* = p*(w) corresponds to a boundary point of F defined by ©. Different boundary 
points can be achieved by choosing different weight vectors in (|20l) . In particular, Lemma \T\ 
implies that q G dF if and only if q = q*(w) for some w > 0. For the analysis in this section, 
it is important to recall from Observation \T\ that there is a bijective map from F onto P and this 
map can be shown to be lfT5l Section 5.3] 



p(q) = (I - G(q)rV)- 1 G(q)rz, q G F 



(23) 



where G(q) := diag(g(gi), . . . ,g(qi()) with g( x ) defined by |(A.3)| and p(G(q)TV) < 1, which 
ensures the existence of p(q) and is satisfied for every q e F liT5l Section 5.3]. Note that g is 
strictly increasing and g(x) > for all x £ Q, which follows from |(A.2~) and |(A.3) 

Lemma 5: q G dF if and only if max ne ^A n (q) = 1 where A n (q) := p(G(q)B^) > 
p(G(q)rV) > 0. 

Proof: The proof is deferred to Appendix IVIII-El ■ 
Now we are in a position to prove the following. 

(g'(qi)/g(qi):- 



Theorem 3: Suppose that |(A.2)] and |(A.3)] hold. Let q G dF and u(q) 
0. Then, we have q = q*(w) given by (1221) if 



g'{qK)/g{q K )) > 



w 



c ■ u(q) o y o x, c > 



(24) 
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where y and x are positive left and right eigenvectors of G(q)B (no) , respectively, associated 
with A no (q) for any n = argmax ngN A n (q). 

Proof: The reader can find the proof in Appendix IVM-Fl ■ 

q-2 = 0(SIR 2 (p)/ 72 ) 




0.5 1.0 1.5 2.0 



a = *(snti(p)/7i) 

Fig. 2. An illustration of Corollary Q] The figure shows an example of the feasible QoS region defined by ifTJ with 2 users 
subject to individual power constraints (C = I). The point q = q*(w) with w = cy o x, c > 0, corresponds to the unique 
max-min SIR-balanced power allocation. The weight vector w is normal to a hyperplane which supports the feasible QoS region 
at q G OF. 

Now we can establish a connection between © and (TfTT) . The connection is illustrated by 
Figure |2] 

Corollary 1: Let y and x be positive left and right eigenvectors of associated with 

p(B (no) ) for any n Q e N (p). If w = cy o x, c> 0, then p = p*(w). 

Proof: The proof can be found in Appendix IVIII-G[ ■ 

B. Saddle point characterization of the Perron roots 

Finally, we point out that Theorem \T\ gives rise to a saddle point characterization of piB^), n 6 
N"o(p). Let U K := {x G R£ : ||x||i = 1} and n+ = Il K n E^ + . We define G : xP + ^l 
as 

C(w,p) := V w k tfj(^ k /SIR k (p)) 
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where ip{x) := —<f)(l/x). A key ingredient in the proof of the saddle point characterization is 
the following theorem, which can be deduced from |[T5l Sections 1.2.4-5]: 

Theorem 4: Assume that |(A. 1 j] - |(A.3)| hold and V is irreducible. Let B = for any n G K, 
and let w = y o x G 11^ where y and x are positive left and right eigenvectors of B associated 
with p(B). Then, for all p > 0, 

*MB))<£^<^). («) 

Equality holds if and only if p = x > 0. 

Now we are in a position to present the saddle point characterization of the Perron root of 
B^, n G N (p), similar to the one for the noiseless case which can be found in lfT5l . 



Theorem 5: Suppose that |(A.2)| and |(A.3)| hold, and V is irreducible. Then, 

^(p(B (no) ))= sup inf G(w,p)= inf sup G(w, p) 
wen+ peP + peP + wen+ 

where n = argmax ne 3sfp(B ( - n ^). Moreover, a point (w*,p*) is the saddle point of G(w,p) if 
and only if p* = p and w G W, where 

W = {w G 11+ : w = c - wW ' ^ c n = 1, c n > 0} 

and := yW o x ( n ) G U+ with p(B {n) )^ = B (n) x("), p(B (n) )y( n ) = (B (n) ) T y(™), 

(y(")) T x (") = l. in words, at the saddle point the power vector is equal to the max-min SIR- 
balanced power allocation, whereas the weight vector is any linear combination of the vectors 

w (n) f()r n e K (p). 

With Theorems \T\ and |4] in hand, the proof is similar to that in lfT5l Section 1.2.4]. The 
existence of a saddle point is ensured by irreducibility of the gain matrix since then positive 
left and right eigenvectors exist. The uniqueness follows from the irreducibility property and the 
normalizations. 

The reason why the theorem is of interest is that it provides a basis for the design of alternative 
power control algorithms for saddle point problems that converge to p and may be amenable to 
distributed implementation. Basically, the idea of the algorithm is redolent of that of primal-dual 
algorithms that converge to a saddle point of the associated Lagrange function ll34l . Development 
of new algorithms is currently a subject of our ongoing work; the main idea, however, consists 
in minimizing the function G(w, p) with respect to p, and simultaneously maximizing G(w, p) 
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with respect to w. The straight-forward approach employs the gradient projection method, where 
in order for the objective function to be convex in the power variable the substitution s = log(p) 
is used. Each iteration encompasses the calculation of the gradient of G(w, e s ), and an update of 
the vectors w and s in the direction and against the direction of the gradient, respectively. This 
process requires a suitable step size to be chosen. The iteration is concluded with a projection 
of the updated values of w and s onto the corresponding sets of valid values. 

The minimization of G(w, p) in the power domain using the gradient projection method cor- 
responds to employing the power control algorithm presented in fTTSl Section 6.5]. In particular, 
the gradient can be computed in a distributed way, and the projection onto the feasible set is also 
distributedly implementable in many cases of interest. As for the optimization in the weights 
domain, the gradient can be computed independently by each node, but performing the projection 
requires in general centralized operation. 

VI. Conclusions 

In this paper, we have characterized the max-min SIR-balanced power allocation under general 
power constraints. This characterization is an extension of the results known previously for 
the noiseless case, in which the power constraints play no role. We have also established a 
connection between the max-min SIR-balancing power allocation problem and the utility-based 
power allocation problem for the considered case, and as an application of our results we have 
discussed two classes of power allocation algorithms based on those two approaches. Finally, we 
have presented a saddle-point characterization of the Perron roots of the extended gain matrices 
which may constitute a basis for developing distributed power control algorithms. 
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VIII. Appendix 

A. Proof of Observation [2] 

Let |(A.3)1 be satisfied, and let p G P+ be any solution to © or, equivalently, ©. Let 



0(SIRfc(p)/7fc), k G X . By Observation [Q F is a convex downward comprehensive set and, 
by |(A.2)[ © and ©, q = (g l7 . . . , q K ) G <9F is its boundary point since at least one power 
constraint is active in the optimum (see Lemma [2]). Thus, by irreducibility of V and Lemma 
[Q there exists w > such that w T (q — u) > for all u G F. Due to positivity of w, this 
implies that for any u G F, u ^ q, there is at least one index i = i(q, u) G % such that 
q-i > Wj. In particular, for any ueF,u/q' there is i = z(q', u) G % such that q[ > where 
q' k = 0(SIRfc(p')/7fc), k G %. On the other hand, however, we have q' < q. This is simply 
because q[ — ■ ■ ■ — q' K — minfc 6 x0(SIRfc(p)). Combining both inequalities shows that q = q', 
and hence, by bijectivity, we obtain p = p', which is unique by Q21 Theorem A.51]. 

B. Proof of Lemma |2] 

Part (i) should be obvious since if we had g n (p) < 1 for all n G N, then it would be 
possible to increase min fee3 <;SIRfc(p)/7fc by allocating the power vector cp G P+ with c = 



1/ ma.x n£ js! g n (p) > 1. In order to show part (ii), note that if |(A.2)| and |(A.3)] hold, then, by 
Observation [Q F is a convex downward comprehensive set. Moreover, p G P+ given by © 
corresponds to a point q G <9F, with % = 0(SIRfc(p)/7fc), k G X. Thus, by irreducibility of V, 
it follows from Lemma \T\ that q is a maximal point of F, and hence q < q for any q G F implies 
that q = q ll34l . That is, there is no vector in F that is larger in all components than q. On 
the other hand, by the discussion in Sect. [Till q is a point where the hyperplane in the direction 
of the vector (1/ K, . . . ,1/ K) intersects the boundary of F. As a result, q x = ■ ■ ■ — q K , which 
together with the maximality property and strict monotonicity of (ft, shows that SIRfc(p)/7fc = (3 
for each k G % where (3 is positive due to (i). If V is irreducible, the uniqueness of p follows 
from Observation |2l 

C. Proof of Lemma |4] 

Let n G N be arbitrary. First we prove part (i). Since 1/P n zc^ > and V is irreducible, 
we can conclude from (fT6l) that B^"- 1 > is irreducible as well. Thus, by the Perron-Frobenius 
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theorem for irreducible matrices (331, Il32lh there exists a positive vector p which is an eigenvector 
of B^ associated with p(B^), and there are no nonnegative eigenvectors of associated 
with p(B^) other than p and its positive multiples. Among all the positive eigenvectors, there 
is exactly one eigenvector p > such that g n (p) = c x . This proves part (i). In order to prove 
(ii), note that if was irreducible, then we could invoke the Perron-Frobenius theorem and 
proceed essentially as in part (i) to conclude (ii) (with the uniqueness property resulting from 
the normalization of the eigenvector so that its last component is equal to c 2 > 0). In order to 
show that A*™- 1 is irreducible, let S(A (n - ) ) be the associated directed graph of {1, . . . , K + 1} 
nodes 11331 . Since TV is irreducible, it follows that the subgraph S(TV) is strongly connected 
11331 . Furthermore, as the vector Tz is positive, we can conclude from (fl4)) that there is a directed 
edge leading from node K + 1 to each node n < K + 1 belonging to the subgraph S(TV). 
Finally, note that as TV is irreducible, each row of TV has at least one positive entry. Hence, 
the vector l/P n c^TV has at least one positive entry as well, from which and (fl4)) it follows that 
there is a directed edge leading from a node belonging to S(rV) to node K + 1. So, S(A^) 
is strongly connected, and thus A^ is irreducible. 



D. Proof of Theorem [7] 

(i) — > (ii): By Lemma [3], p G P + satisfies (fT5l) for some (3 > 0. Thus, by Lemma H part (i) 
implies part (ii). (ii) — -> (Hi): Given any n G No(p)> it follows from (fT5l) that p(B ( - n ^)p = B^p 
with g n (p) = 1 is equivalent to p(B^)p = TVp + Tz, which in turn can be rewritten to give 
(fT3l with (3 = pCB^) and p = (p, 1). Since p is positive, so is also p. Thus, p with Pk+i = 1 
is a positive right eigenvector of A*-™- 1 and the associated eigenvalue is equal to p(B^) > 0. So, 
considering part (Hi) of Lemma H we can conclude that (Hi) follows from (ii). (Hi) — > (i): By 
Lemma HI for each n e No(p)> there exists exactly one positive vector p with px+i = 1 such 
that (IT~3b is satisfied. Furthermore, j3 = p(A ( - n ^ > ), n £ 3Sf (p) is a simple eigenvalue of A^ n \ Now 
considering Lemma [3] proves the last missing implication. 



E. Proof of Lemma \5\ 



By © with |(A.2)| we have q G F if and only if there is p G P such that 0(SIRfc(p)/7fc) > qu 
for each k G X. Thus, q G F if and only if 1/A := max peP mm ke x(^^Rk(p)/lk9(lk)) > 1 
where the maximum always exists. Comparing the left hand side of the inequality above with 
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© shows that the only difference to the original problem formulation is that 7 fc is substituted by 
Ikdilk) or, equivalently, T by G(q)T, which is positive definite as well. Thus, by (flOl) . Theorem 
[T]and Theorem [2} we have q G F if and only if A = max ne:N - A n (q) < 1. Moreover, p(q) given 
by (1231 is the unique power vector such that q k = 0(SIRfc(p(q))/7/,.) for each k G %. Since the 
Neumann series converges for any q G F, we have p(q) = ^^ (G(q)rV) / G(q)rz. Now as 
G(q)Tz is positive and G(q)TV is irreducible, we can conclude from |fT51 Lemma A. 28] and 



(A.2)| that each entry of p(q) is strictly increasing in each entry of q. Thus, as F is downward 



comprehensive and q ^ <9F holds if and only if all power constraints are inactive, for every 
q G int(F), there is q G <9F such that q = q + u for some u > 0. By irreducibility of B^, this 
implies that A n (q) < A n (q + u) = A n (q) < 1 for each n G N So, if max ne jj- A n (q) = 1, then 
q G <9F. Conversely, if q G <9F, we must have max„ e ^A n (q) = 1 since otherwise there would 
exist q ^ F such that max ng ^A n (q) = 1, which would contradict Theorem [2l This completes 
the proof. 



F. Proof of Theorem \3\ 

Let q G dF and n = argmax„ g ^ A n (q) be arbitrary and note that F is a convex set. So, by 
Lemma [U there is w > such that q maximizes q i— > w T q over F. Lemma [5] implies that this 
convex problem can be stated as q = argmax q w T q subject to A no (q) = 1, q G Q K . Due to 
(A.2)[ the spectral radius is continuously differentiable on Q K . Thus, the Karush-Kuhn-Tucker 
conditions ll34l . which are necessary and sufficient for optimality (due to the convexity property), 
imply that w is parallel with VA no (q). Now, by [42], we have = y k g'(q k ) Eiex^^i = 

g ^ykI2iex9(<lk)b k J 'xi = X no (q)^^y k x k = ^jVkXk for each k G % where y and x are 
left and right positive eigenvectors of G(q)B^ n °' 1 associated with A no (q), which, by irreducibility, 
are unique up to positive multiples. 



G. Proof of Corrolary \T\ 

As V is irreducible, Observations [T] and [2] (see also the proof) imply that p corresponds to a 
point q G dF. Since q*(w) G dF for any w > 0, it follows from Theorem [3] that q = q*(w) if 
w is has the form (l24l) . Now by Observations [TJ and [2l we have qi — • ■ ■ — q~K- Thus, as both g 
and its derivative g' are strictly monotonic (by |(A.2)1 and |(A.3")] ), we must have u(q) = al, a > 
and G(q) = l/p(B ( - no ^)I. Thus, the corollary follows from Theorem [3] 
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